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SUMMARY 

Random  sea  surfaces  representative  of  sea  states  3 through  8 have  been 
generated  for  use  with  the  computer  program  CABUOY.  Output  waveforms  are  pre- 
sented graphically,  along  with  the  frequencies,  eunplitudes  and  phases  of  the 
component  waves.  Physical  characteristics  of  the  generated  waveforms  are  com- 
pared with  actual  sea  state  data. 
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I BACKGROUND 

f C.vBUOY  is  a Fortran  computer  program  which  solves  for  the  two  dimensional 

i dynamic  motions  of  ocean  deployed  cable  systems.^  The  program  calculates  the 

I response  of  a moored,  drifting  or  towed  system  to  a surface  wave  forcing  func- 

i tion,  and  outputs  the  position  of  each  suspended  body  and  the  resultant  cable 

I tensions.  Designed  as  a developmental  tool  for  the  evaluation  of  new  sonobuoy 

systems,  CABUOY  requires  a standardized  input  for  surface  motions  so  that  en- 
gineering changes  can  be  compared  against  known  designs.  Sea  surfaces  repre- 
sentative of  sea  states  3 through  8 were  desired. 
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DISCUSSION 

CABUOY  generates  a random  surface  motion  through  the  superposition  of 
sine  waves  of  the  form 

y = A sin  (kx  - ut  + 6) 

where 

y = vertical  coordinate  of  wave 
A = wave  amplitude 
k = 2it/X  = wave  number 
X = wave  length 

X = horizontal  coordinate  of  wave 
0)  = 2Tr/T  = circular  frequency 
T = wave  period 
t = time 
0 = wave  phase 

This  effort  was  concerned  with  the  selection  of  component  amplitudes  and 
periods  which  would  yield  motions  characteristic  of  wind  driven  waves  at  sea, 
in  terms  of  significant  wave  height  and  period.  As  CABUOY  is  limited  to  19 
component  waves,  the  pattern  would  not  be  completely  random.  It  was  considered 
that  1000  seconds  of  random  motion  would  be  sufficient  in  most  cases. 

The  calculation  of  component  wave  amplitudes  was  based  upon  a knowledge 
of  the  energy  content  of  wind  driven  waves.  According  to  the  Pierson-Moskowitz 
wave  theory,^  the  energy  spectral  density  for  a fully  developed  sea  can  be  de- 
fined by 


S(a))  = 8.1  X 10 


■3  2 -5 

g to 


^-.74(g/M^to) 


1.  A Fortran  IV  Computer  Program  for  the  Time  Domain  Analysis  of  the  Two 
Dimensional  Dynamic  Motions  of  General  Buoy-Cable-Body  Systems , H.  T. 
Many,  DTNSRDC  Report  77-0046,  June  77. 

2.  Buoy  Engineering , H.  O.  Berteaux,  John  Wiley  & Sons,  1976. 
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where 

g = acceleration  due  to  gravity 
u = wind  velocity 

Solutions  to  this  equation  are  illustrated  in  figure  1 for  several  wind 
velocities.  Integration  of  this  curve  for  a constant  wind  velocity  yields  the 
cocumulative  energy  function,  illustrated  in  figure  2.  This  function  repre- 
sents the  energy  content  of  a fully  developed  sea  as  a function  of  wave  period. 
For  a constant  wind  velocity,  the  energy  defined  by  the  cocumulative  function 
at  a period  To  represents  the  total  energy  for  all  waves  of  periods  To  or  less. 
Therefore,  for  some  interval  AT,  the  quantity 

E(To  + AT/2)  - F.(To  - AT/2) 

represents  the  total  wave  energy  contained  in  a period  band  centered  at  To  with 
a width  of  AT.  Thus,  from  the  cocumulative  spectrum,  it  is  possible  to  calcu- 
late the  total  wave  energy  existing  in  discrete  period  bands. 

From  statistical  theory,  the  following  relationships  between  wave  height 
and  wave  energy  have  been  defined:^ 

Significant  wave  height  = 2 . 83'\/energy 
Tenth  highest  wave  = 3. 60^ energy 

Average  wave  height  = 1 .77'\/ energy 

For  this  analysis,  the  component  wave  heights  were  defined  according  to 
the  definition  for  average  wave  height  using  the  energy  derived  from  the  co- 
cumulative spectrum.  It  was  found  that  component  waves  of  average  height 
summed  together  generated  sea  states  consistent  with  recorded  data,  as  summa- 
rized in  the  sea  state  chart,  table  I, 

Component  frequencies  were  determined  by  dividing  the  significant  period 
band,  defined  as  the  band  of  periods  on  the  cocumulative  distribution  function 
that  excludes  the  top  5 percent  and  lower  3 percent  of  the  total  energy,  into 
equally  spaced  intervals.  The  component  waves  were  then  assigned  the  period 
in  the  center  of  the  interval.  Component  periods  were  shifted  slightly,  if 
necessary,  to  avoid  harmonics  which  would  decrease  the  randomness  of  the  re- 
sultant wave  pattern. 

Tables  II  through  VII  list  the  component  wave  frequencies  and  amplitudes 
as  calculated  for  sea  states  3 through  8.  The  phase  shifts  were  included  to 
eliminate  the  unnaturally  large  wave  which  resulted  when  all  components  were 
summed  together  beginning  at  time  0.  The  phase  shift  amounts  to  shifting  the 
initial  time  20  seconds,  and  eliminates  the  singularity.  One  thousand  seconds 
of  random  sea  state  data  can  be  generated  before  the  pattern  repeats. 


3,  wind  Waves,  Their  Generation  and  Propagation  on  the  Ocean  Surface, 
B.  Kinsman,  Prentice  Hall,  1965. 
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CO-CUMULATIVE  ENERGY  FUNCTION 
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FIGURE  2 - Co -Cumulative  Energy  Function 
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T A B 1,  I;  II 

SINK  WAVE  COMl’ONENT  INPUTS  FOR  SEA  STATE  GENERATOR 
SEA  STATE  3 


Frequency 

Amplitude 

Phase 

(Hz) 

(ft) 

(deg) 

0.459 

0.090 

64.8 

0.394 

0.127 

316.8 

0.345 

0.270 

324.0 

0.307 

0.238 

50.4 

0.277 

0.324 

194.4 

0.252 

0.382 

14.4 

0.231 

0.312 

223.2 

0.213 

0.402 

93.6 

0.198 

0.349 

345.6 

0.185 

0.349 

252.0 

0.174 

0.493 

172.8 

0.164 

0.493 

100.8 

0.154 

0.285 

28.8 

0.146 

0.285 

331.2 

0.139 

0.201 

280.8 

0.133 

0.201 

237.6 

0.127 

0.285 

194.4 

0.121 

0.201 

151.2 

0.116 

0.201 

115.2 
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TABLE  III 

SINE  WAVE  COMPONENT  INPUTS  FOR  SEA  STATE  GENERATOR 
SEA  STATE  4 


Frequency 

Amplitude 

Phase 

(Hz) 

(ft) 

(deg) 

0.371 

0.255 

151.2 

0.323 

0.238 

165.6 

0.287 

0.349 

266.4 

0.258 

0.402 

57.6 

0.234 

0.493 

244.8 

0.214 

0.349 

100.8 

0.197 

0.604 

338.4 

0.183 

0.569 

237.6 

0.171 

0.402 

151.2 

0.160 

0.569 

72,0 

0.151 

0.636 

7.2 

0.142 

0.402 

302.4 

0.135 

0.569 

252.0 

0.128 

0.402 

201.6 

0.122 

0.285 

158.4 

0.116 

0.285 

115.2 

0.111 

0.285 

79.2 

0.106 

0.569 

43.2 

0.069 

0.402 

136.8 
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TABLE  IV 

SINE  WAVE  COMPONENT  INPUTS  FOR  SEA  STATE  GENERATOR 
SEA  STATE  5 


Frequency 

Amplitude 

Phase 

(Hz) 

(ft) 

(deg) 

0.253 

0.604 

21.6 

0.224 

0.667 

172.8 

0.2 

0.753 

0.0 

0.182 

0.944 

230.4 

0.166 

0.9 

115.2 

0.153 

1.026 

21.6 

0.142 

1.065 

302.4 

0.132 

0.986 

230.4 

0.124 

1.365 

172.8 

0.116 

1.102 

115.2 

0.110 

0.805 

36.0 

0.104 

1.026 

28.8 

0.099 

0.854 

352.8 

0.094 

0.9 

316.8 

0.089 

0.697 

280.8 

0.086 

0.753 

259.2 

0.082 

0.753 

230.4 

0.079 

0.493 

208.8 

0.076 

0.402 

187.2 
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TABLE  V 

SINE  WAVE  COMPONENT  INPUTS  FOR  SEA  STATE  GENERATOR 
SEA  STATE  6 


Frequency 

Amplitude 

Phase 

(Hz) 

(ft) 

(deg) 

0.209 

0.944 

64.8 

0.186 

1.138 

259.2 

0.168 

1.102 

129.6 

0.153 

1.273 

21.6 

0.141 

1.273 

295.2 

0.130 

1.365 

216.0 

0.121 

1.777 

151.2 

0.113 

1.635 

93.6 

0.106 

1.559 

43.2 

0.1 

1.304 

0.0 

0.095 

1.61 

324.0 

0.09 

1.479 

288.0 

0.085 

0.9 

252.0 

0.081 

1.207 

223.2 

0.078 

1.173 

201.6 

0.074 

1.102 

172.8 

0.071 

0.9 

151.2 

0.068 

0.9 

129.6 

0.066 

1.273 

115.2 
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TABLE  VI 

SINE  WAVE  COMPONENT  INPUTS  FOR  SEA  STATE  GENERATOR 
SEA  STATE  7 


Frequency 

Amplitude 

Phase 

(Hz) 

(ft) 

(deg) 

0.171 

1.423 

151.2 

0.153 

1.635 

21.6 

0.139 

1.660 

280.8 

0.127 

2.347 

194.4 

0.117 

2.205 

122.4 

0.108 

2.546 

57.6 

0.101 

2.381 

7.2 

0.094 

2.546 

316.8 

0.088 

2.514 

273.6 

0.083 

2.432 

237.6 

0.079 

2.546 

208.8 

0.075 

2.530 

180.0 

0.071 

2.012 

151.2 

0.068 

2.012 

129.6 

0.065 

1.684 

108.0 

0.062 

1.909 

86.4 

0.060 

1.273 

72.0 

0.057 

1.559 

50.4 

0.055 

1.273 

36.0 
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TABLE  VII 

SINE  WAVE  COMPONENT  INPUTS  FOR  SEA  STATE 
SEA  STATi:  8 


Frequency 

Amplitude 

(Hz) 

(ft) 

0.135 

2.347 

0.121 

3.509 

0.110 

3.745 

0.101 

3.656 

0.093 

3.84 

0.086 

3.711 

0.080 

4.316 

0.075 

4.5 

0.071 

4.124 

0.067 

4.316 

0.063 

4.589 

0.060 

3.711 

0.057 

4.124 

0.054 

3.486 

0.052 

3.367 

0.050 

2.846 

0.048 

2.7 

0.046 

1.8 

0.044 

1.559 

GENERATOR 


Phase 

(deg) 

252.0 

151.2 
72.0 

7.2 

216.0 

259.2 
216.0 
180.0 

151.2 
122.4 

93.6 

72.0 

50.4 
28.8 

14.4 

0.0 

345.6 

331.2 
316.8 
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Figures  3 through  8 present  600  seconds  of  wave  data  for  sea  states  3 
through  8 which  were  generated  using  the  previously  defined  component  waves. 
These  waves  coinpare  favorably  with  the  sea  state  chart  in  table  I.  Figure  9 
presents  actual  wave  data  recorded  at  sea  north  of  St.  Croix  in  the  U.S.  Virgin 
Islands.  In  general,  the  confuted  sea  data  compare  favorably  with  actual  data. 


CONCLUSION 

Random  sea  surfaces  composed  of  discrete  sine  waves  whose  amplitudes  and 
frequencies  are  derived  from  the  cocumulative  sea  spectrum  can  be  generated 
which  compare  favorably  to  actual  sea  surface  date.  These  computer  sea  states 
should  be  used  with  the  program  CABUOY  when  it  is  necessary  to  evaluate  the  re- 
sponse of  a cable  suspended  system  to  a random  sea  state  forcing  function. 


FIGURE  3 - Computer  Generated  Sea  State  3 Wave 
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STATE  DATA 


FIGURE  9 - Measured  Sea  State  3 Data 
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